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ABSTRACT 


An  experimental  and  theoretical  investigation  was  conducted  to 
determine  the  degree  of  jet  boundary  simulation  obtainable  for  under¬ 
expanded  jets  exhausting  into  a  quiescent  atmosphere.  Tests  were 
conducted  with  nitrogen,  carbon  dioxide,  and  helium  gases,  and 
theoretical  boundaries  were  obtained  by  a  method- of- characteristics 
solution.  It  was  determined  that  the  method- of- characteristics  solu¬ 
tion  accurately  represents  the  experimental  jet  boundaries  and  that 
the  matching  of  the  parameters  6j  and  Mj/y  gives  good  boundary 
simulation  for  a  wide  range  of  jet  conditions.  Also  a  rapid,  simple 
method  for  estimating  jet  boundary  shape  is  developed  from  the 
method-of- characteristics  results. 
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SECTION  I 
INTRODUCTION 


In  recent  years,  much  interest  has  arisen  in  the  structure  and 
properties  of  highly  underexpanded  rocket  or  ramjet  exhaust  plumes. 
This  interest  has  been  created  by  problems  such  as  heating  and  eros¬ 
ion  of  adjacent  surfaces,  radar  and  communication  signal  interference, 
separation  of  flow  over  the  vehicle,  and  radiation  emitted  by  the  large 
exhaust  plumes.  Studies  of  and  attempts  to  obtain  solutions  to  these 
problems  have  led  to  the  construction  of  high  altitude  test  facilities 
and  the  adaptation  of  some  existing  facilities.  Methods  used  for  the 
experimental  simulation  of  the  full-scale  jets  vary  in  complexity  from 
the  use  of  small  cold-gas  jets  to  an  almost  exact  duplication  of  the  full- 
scale  jets.  The  degree  of  similitude  required  and/or  used  depends  on 
the  particular  problem  under  investigation  and  is  discussed  in  Ref.  1. 
However  the  tremendous  vacuum  pumping  capacities  required  and  the 
complexity  and  cost  of  hot  rocket  tests  frequently  dictate  that  studies 
be  conducted  with  small  jets  of  cold  gases. 

One  jet  characteristic  frequently  desired  is  that  the  boundary  of 
the  experimental  jet  duplicates  that  of  the  full-scale  jet.  In  an  early 
study.  Love  (Ref.  2)  concluded  that  for  slightly  underexpanded  jets 
(Pj/  Poo  from  1  to  10)  good  boundary  simulation  was  obtained  if  the 
initial  angle,  6 j,  of  the  jet  boundary  was  matched.  Goethert  (Ref.  3) 
later  proposed  matching  the  parameters  5j  and  -yM for  highly  under¬ 
expanded  jets.  Pindzola  (Ref.  4)  concluded  that  better  simulation  could 
be  obtained  by  matching  6 j  and  yM \^l$i  where  the  Mach  number  is  now 
based  on  the  jet  boundary  conditions  rather  than  the  nozzle  exit  condi¬ 
tions  as  proposed  by  Goethert. 

The  objective  of  the  investigation  reported  herein  was  to  determine 
the  degree  of  jet  boundary  simulation  obtained  using  Pindzola's  param¬ 
eters  and  to  determine  if  other  simulation  parameters  might  prove  more 
useful.  The  investigation  included  both  experimental  and  theoretical 
studies  of  the  jet  boundaries.  Values  of  the  simulation  parameters  6j 
and  were  selected  that  were  representative  of  underexpanded 

jet  plumes,  and  tests  were  conducted  with  nitrogen  (N2),  helium  (He), 
and  carbon  dioxide  '(CO2)  jets  at  constant  values  of  the  parameters. 

In  addition,  jet  boundaries  were  calculated  for  several  sets  of  simu¬ 
lation  parameters  by  an  approximate  method  and  also  by  a  method-of- 
characteristics  solution. 

From  these  studies  it  was  determined  that  although  the  parameters 
6j  and  yMi^l  fix  gave  good  jet  boundary  duplication,  better  duplication 
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could  be  obtained  by  matching  the  parameters  6j  and  Mi/ 7.  Also,  a 
simple  method  was  developed  for  rapidly  predicting  the  jet  boundary 
for  a  wide  range  of  conditions. 


SECTION  II 

THEORETICAL  B0LT7ARY  CALCULATIONS 


2.1  LATVALA'S  APPROXIMATION 

The  approximation  used  by  Pindzola  (Ref.  4)  for  the  original 
evaluation  of  the  simulation  parameters  6j  and  yM^lfii  was  devel¬ 
oped  by  Latvala  in  Ref.  5.  Latvala's  approximation  is  based  on  quasi- 
one -dimensional  relations  and  the  assumption  of  isentropic,  radial 
flow  from  the  nozzle.  It  is  an  adaptation  of  a  method  presented  by 
Adamson  and  Nichols  in  Ref.  6,  but  uses  a  spherical  rather  than  a 
planar  area  to  define  the  flow  area  of  the  jet. 


The  approximate  solution  starts  with  the  known  initial  angle  of 
the  jet,  6 j,  calculated  from 


8 


j 


=  v 


1 


e 


N 


(1) 


The  condition  of  constant  pressure  along  the  jet  boundary  is  then 
used  to  determine  the  changes  in  jet  boundary  angle  required  to  com¬ 
press  the  flow  and  balance  the  pressure  decrease  caused  by  the  one¬ 
dimensional  flow  area  increase.  The  relation  between  turning  angle 
and  pressure  change  is  given  by  the  Prandtl- Meyer  expansion  relation 


(Ay)  +  higher  order  terms  in 


(2) 


which  is  identical  (except  for  sign  of  Ai/)  to  the  expression  for  a  com¬ 
pression  through  an  oblique  shock  wave,  up  to  and  including  the  term 
in  (Av)2.  It  is  apparent  then  that  two  jets  that  have  the  same  value  erf 
the  coefficient  yMy^lfii  could  be  expected  to  have  the  same  pressure 
change -deflection  angle  relation  within  the  limits  of  the  linearization  of 
Eq.  (2).  This  therefore  gives  rise  to  the  use  of  Sj  and  yM^/fii  as  the 
jet  boundary  simulation  parameters. 


2.2  METHOD-OF-CHARACTERISTICS  SOLUTION 

The  most  accurate  method  existing  for  calculating  the  properties 
of  an  inviscid  expanding  jet  is  the  method  of  characteristics.  However, 
the  accuracy  obtained  in  the  actual  application  depends  on  factors  such 
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as  the  knowledge  of  the  actual  flow  properties  at  the  start  of  the  cal¬ 
culation  (the  nozzle  exit),  the  frequent  assumption  of  neglecting  the 
jet  boundary  shock,  and  the  nonideal  behavior  of  the  gas.  Also,  the 
accuracy  is  limited  by  the  computer  program  used  because  of  the  finite 
mesh  size  rapidly  becoming  larger  as  the  solution  proceeds  downstream. 
In  most  applications,  the  effects  of  the  viscous  mixing  on  the  jet  bound¬ 
ary  must  be  considered. 

Application  of  method -of- characteristics  solutions  to  determine 
the  free  jet  boundaries  have  been  made  by  several  authors  (Refs.  2,  5, 
and  7  through  12).  Among  these,  Refs.  2,  5,  and  7  included  limited 
experimental  data  for  sonic  nozzles  or  low  to  moderately  expanding 
jets  from  supersonic  nozzles.  Experimental  determination  of  jet 
boundaries  for  highly  underexpanded  jets  has  been  prevented  by  the 
density  limitations  of  conventional  optical  instrumentation;  therefore 
a  "glow  discharge"  technique  was  used  in  Ref.  12.  The  method-of- 
characteristics  solution  was  also  compared  in  Ref.  13  with  experi¬ 
mental  data  on  the  internal  structure  of  jets  expanding  into  a  vacuum. 

In  general,  the  above  authors  found  the  method  of  characteristics  to 
reasonably  predict  the  jet  boundaries  and  flow  properties,  within  its 
inherent  limitations.  A  good  survey  of  some  of  the  exact  and  approxi¬ 
mate  solutions  for  the  jet  structure  has  been  given  by  Adamson 
(Ref.  14). 

The  computer  program  used  for  the  present  calculations  was 
developed  by  the  Lockheed  Missile  and  Space  Corporation  for  use  with 
the  IBM  7090  computer.  1  This  program  computes  the  supersonic  flow 
properties  of  a  gas  discharging  from  an  underexpanded  nozzle  into  either 
a  quiescent  atmosphere  or  a  hypersonic  free  stream.  Ideal  gas  relations, 
oblique  shock  relations  for  the  boundary  shock,  and  the  axisymmetric 
method  of  characteristics  are  employed.  A  complete  description  of  the 
program  is  given  in  Ref.  15,  and  only  the  principal  assumptions  and  the 
techniques  used  in  the  application  of  the  program  are  given  here. 

Jet  boundaries  presented  in  this  report  were  computed  on  an  IBM 
360/50  computer  for  the  following  conditions  and  assumptions: 

1.  For  the  case  of  a  quiescent  atmosphere  only. 

2.  Nozzle  exit  flow  properties  were  computed  by  the 
program  from  ideal  gas,  one -dimensional  flow 
relations  utilizing  point  source  flow  (spherical  areas). 

Specific  heat  ratios  of  1.  28,  1.  40,  and  1.  667  were 
used  for  both  the  nozzle  and  the  jet  flows.  The  value 
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of  1.  28  used  for  comparison  with  experimental  tests  of 
CO2  corresponds  to  excitation  of  one  of  the  four  modes 
of  the  molecular  vibration. 

3.  The  method- of- characteristics  computation  starts  on  the 
spherical  nozzle  exit  surface;  normally  4Q,  points  were 
used  for  the  origin  of  the  characteristic  net. 

4.  FLow  properties  for  the  initial  expansion  around  the  nozzle 
lip  are  computed  from  Prandtl-Meyer  relations.  Approxi¬ 
mately  1.  5-deg  increments  were  used  for  the  Prandtl- 
Meyer  expansion  system. 

Recently  a  new  program  has  been  developed  (Ref.  16)  which  in¬ 
cludes  equilibrium  reacting  gas  mixtures  and  a  wider  variety  of  flow 
field  problems.  A  comparison  of  this  program  with  some  experimental 
data  is  given  in  Ref.  17.  These  comparisons  include  limited  data  from 
jets  with  sonic  nozzles  and  with  low  exit  pressure  ratio  supersonic 
nozzles.  In  general  the  comparisons  were  considered  very  good.  This 
program  (Ref.  16)  is  currently  being  adapted  to  the  AEDC  computer 
system. 

2.3  JET  BOUNDARY  VISCOUS  MIXING  REGION 

The  velocity  gradients  normal  to  the  flow  direction  in  the  interior 
of  the  jet  are  usually  small,  and  therefore  viscosity  effects  are  neg¬ 
lected.  The  jet  interior  structure  can  then  be  calculated  by  the  method 
of  characteristics  as  described  previously.  However,  the  gradients  are 
not  negligible  at  the  jet  boundary,  and  effects  of  viscous  mixing  must  be 
considered. 

Ideally,  the  inviscid  boundary  can  be  determined  and  the  flow  proper¬ 
ties  at  the  inviscid  boundary  used  to  determine  the  corresponding  mixing 
region  flow  properties.  The  viscid  mixing  solution  can  then  be  super¬ 
imposed  on  the  inviscid  jet  solution  to  define  the  physical  jet.  This  pro¬ 
cedure  is  probably  valid  if  the  mixing  region  is  small  compared  to  the 
jet  radius  and  if  no  chemical  reactions  occur. 

The  growth  of  the  turbulent  mixing  region  was  calculated  for  the 
conditions  corresponding  to  the  present  experimental  tests  using  the 
method  given  by  Bauer  in  Ref.  18.  The  type  of  mixing  considered  is 
turbulent,  two-dimensional,  compressible,  isobaric,  isoenergetic, 
and  without  an  initial  boundary  layer.  It  is  of  significance  to  note  that 
for  the  high  jet  boundary  Mach  numbers  considered  in  this  report,  the 
mixing  region  thickness  and  velocity  distributions  are  almost  identical 
for  all  jets.  For  these,  the  edges  of  the  mixing  regions  are: 
n/re  /&i/re  =  0. 13  for  the  outer  edge  and  -0.  01  for  the  inner  edge. 
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SECTION  III 

APPARATUS  AND  TEST  PROCEDURE 

3.1  TEST  HARDWARE  AND  PROCEDURE 


The  experiments  reported  in  this  report  were  conducted  in  an 
18-in.  test  cell  connected  to  a  six-stage  steam  ejector.  The  test  cell 
installation  is  shown  in  Fig.  1,  Appendix  I.  Test  gases  were  heated 
by  a  resistance -type  heater  and  expanded  through  conical  nozzles 
into  the  test  cell.  Test  nozzles  had  throat  diameters  of  approximately 
0.  1  in. 


Two  sets  of  values  of  the  simulation  parameters  were  selected 
as  being  representative  of  typical  underexpanded  jet  plumes: 


and 


8 j  =  60  deg,  yMj/jB,  =  14 
=  60  deg,  =  16 


These  values  of  the  simulation  parameters  were  maintained  constant, 
and  tests  were  conducted  with  pure  gases  in  10-  and  20-deg  half-angle 
conical  nozzles.  Nitrogen,  CO2,  and  He  gases  were  tested  at  a  simu¬ 
lation  parameter  value  of  14;  N2  and  He  were  tested  at  a  value  of  18. 
Different  values  of  ambient  pressure  and  nozzle  area  ratio  were 
required  for  each  of  10  different  nozzles  in  order  to  maintain  the 
selected  values  of  the  simulation  parameters. 


Since  each  gas  required  a  different  nozzle  area  ratio  to  match  the 
selected  value  of  the  simulation  parameters,  testing  began  with  CO2 
(the  largest  nozzle  area  ratio).  The  nozzle  exit  was  then  cut  off  to 
give  the  required  area  ratios  for  N2  and  He.  The  configurations  of 
the  various  test  nozzles  are  shown  in  Fig.  2. 


Possible  effects  of  nozzle  boundary  layer,  jet  viscous  mixing, 
gas  condensation,  specific  heat  ratio  variation,  and  gas  nonequilibrium 
could  not  be  precisely  evaluated  beforehand;  therefore,  tests  were 
conducted  over  a  range  of  stagnation  pressures  from  100  to  400  psia;  and, 
temperatures  from  100  to  1000°F.  At  each  test  condition,  cell  ambient 
pressure  was  adjusted  to  give  the  required  jet  expansion  to  match  the 
chosen  values  of  the  simulation  parameters.  All  tests  were  conducted 
under  steady-state  conditions. 
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3.2  FLOW  VISUALIZATION 

A  glow  discharge  technique  was  used  to  illuminate  the  jet  plume 
so  that  it  could  be  photographed,  since  the  required  test  cell  pres¬ 
sures  were  well  below  the  sensitivity  limit  of  conventional  optical 
techniques.  The  glow  discharge  technique,  previously  used  for  jet 
plume  studies  in  Ref.  12,  utilized  an  electrically  charged  probe 
approximately  5  in.  downstream  of  the  nozzle  exit.  This  electrical 
charge  was  supplied  by  Os  ran#  ignition  unit  originally  designed  for  igni¬ 
tion  of  2000-w  high-pressure  xenon  arc  lamps.  This  unit  is  rated  at 
40  kv-RF  at  18  ma. 

Very  small  amounts  of  N2  or  CO2  gas  were  bled  into  the  upstream 
end  of  the  test  cell  through  a  manifold  ring.  This  technique,  discussed 
in  Ref.  12,  provides  a  fine  adjustment  of  the  test  cell  pressure  but  was 
used  primarily  to  give  a  color  contrast  between  the  jet  and  the  ambient 
atmosphere.  This  is  achieved  by  inbleeding  CO2  when  using  N2  as  the 
jet  gas  and  vice  versa.  Although  the  color  of  the  gas  discharge  varies 
somewhat  with  pressure  level  and  gas  temperature,  the  generally  blue 
glow  of  the  CO2  improves  the  boundary  contrast  between  the  pink  glow 
of  the  N2.  Carbon  dioxide  was  also  used  for  the  inbleed  gas  in  the 
tests  with  He.  Figure  3  gives  an  interior  view  of  the  test  cell  showing 
the  electrical  probe  and  the  gas  inbleed  manifold  installation. 

The  illuminated  plumes  were  photographed  with  a  3 5 -mm  camera 
equipped  with  a  50-mmlens.  Ektachrome®  ER-ASA1 60 -Tungsten 
Balance  film  was  used,  and  exposure  times  were  normally  1/30  sec. 

Data  were  obtained  at  several  exposure  settings  at  each  test  condi¬ 
tion,  with  f/4  generally  giving  best  results.  Optimum  exposure  varied 
with  jet  stagnation  temperature  and  pressure.  No  change  in  boundary 
location  could  be  detected  as  a  result  of  a  change  in  exposure. 

Typical  photographs  of  the  jet  plumes  are  shown  in  Fig.  4.  Also 
shown  in  the  figure  are  the  boundary  points  that  were  selected  as  the 
experimental  data  as  well  as  the  theoretical  method-of -characteristics 
solutions.  The  original  color  positives  were  projected,  and  both  top 
and  bottom  boundaries  were  measured  and  averaged  to  give  the  plotted 
data  presented  later  in  this  report. 

The  light  emitted  by  the  excited  gas  has  been  shown  to  be  nearly  a 
linear  function  of  the  gas  density  for  the  pressure  range  of  these  tests 
(Ref.  19).  For  the  highly  expanded  plumes  investigated,  the  gas  density 
very  rapidly  decreases  downstream  of  the  nozzle  exit  as  typically  shown 
in  Fig.  5.  However,  downstream  at  the  boundary  shock  (which  lies  close 
to  the  inviscid  boundary)  an  order -of -magnitude  rise  in  density  occurs. 
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giving  the  bright  band  of  light.  The  outer  edge  of  this  band  was  meas¬ 
ured  as  being  the  jet  boundary.  Near  the  nozzle  exit  the  jet  boundary 
tends  to  disappear,  particularly  for  the  CO2  and  the  N2  jets.  The 
disappearance  of  the  boundary  may  be  attributed  to  a  number  of  effects, 
including  the  shorter  viewing  path  through  the  three-dimensional  jet 
at  the  smaller  axial  distances,  glow  phenomena  electrode  effects,  and 
the  nozzle  boundary  layer. 

3.3  PRESSURE  PROBE  TRAVERSES 

To  confirm  that  the  visible  boundary  did  in  fact  closely  represent 
the  physical  jet,  limited  pressure  probe  traverses  were  made  through 
the  CO2  and  N2  jet  'boundaries.  These  traverses  were  made  parallel 
to  the  jet  axis  and  do  not  necessarily  give  the  pitot  pressure.  The 
probe  tip  was  1/8  in.  in  diameter  and  is  shown  in  Fig.  3.  Results  of 
these  tests  are  shown  in  Fig.  6.  Generally  good  agreement  is  shown 
between  the  jet  boundary  indicated  by  the  pressure  probe  traverse  and 
that  obtained  with  the  glow  discharge  technique. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


4.1  THEORETICAL  JET  BOUNDARIES 

Jet  boundaries  were  calculated  by  the  Latvala  approximation  and 
by.  the  method  of  characteristics  for  values  of  the  simulation  param¬ 
eters  of  =  60  deg  and  equal  to  10,  14,  18,  and  22. 

Boundaries  were  computed  corresponding  to  nozzle  half-angles  of  10 
and  20  deg  with  jet  specific  heat  ratios  of  1.667,  1.40,  and  1.  28. 
These  calculations  cover  a  very  wide  range  oi  possible  operating 
conditions,  as  shown  in  Table  I  (Appendix  ID  and  in  Fig.  7. 

Results  of  the  calculations  by  the  Latvala  approximation  are 
shown  in  Fig.  8,  and  by  the  method  of  characteristics  in  Fig.  9. 

Note  that  since  the  Latvala  approximation  starts  after  the  jet  expan¬ 
sion  to  ambient  pressure,  the  calculations  shown  in  Fig.  8  are  valid 
for  all  values  of  the  nozzle  half -angle,  These  computations  show 

the  Latvala  approximation  to  give  almost  exact  boundary  duplication 
for  all  conditions  calculated,  whereas  the  method  of  characteristics 
predicts  substantial  divergence  between  boundaries,  especially  at  the 
lower  values  of  7Mi2/j3i.  Somewhat  closer  boundary  duplication  for 
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the  method- of- characteristics  solutions  is  obtained  by  adding  the  nozzle 
half-angle,  0^,  as  an  additional  parameter.  Also  note  that  the  bound¬ 
aries  with  a  specific  heat  ratio,  y,  of  1.  667  diverge  the  greatest,  and 
that  for  most  applications  good  boundary  simulation  would  be  obtained 
for  axial  distances  of  10  to  15  nozzle  exit  radii  by  using  air  (y  =  1.4) 
to  simulate  a  hot  rocket  (y  =  1.  28). 


4.2  EXPERIMENTAL  RESULTS 

Results  of  experimental  tests  with  CC>2,  N2,  and  He  for  yM =  14 
are  shown  in  Fig.  10,  and  tests  with  N2  and  He  for  yM  J2/ jSj  =  18  are  shown 
in  Fig.  11.  Carbon  dioxide  was  not  used  in  the  tests  with  =  is 

because  of  the  very  low  ambient  pressures  and  the  very  large  nozzle 
area  ratios  (Fig.  7)  required.  Data  symbols  are  only  shown  for  the 
portions  of  the  boundaries  for  which  a  distinct  boundary  could  be  deter¬ 
mined  from  the  projected  negatives  and  are  shown  as  read,  without 
any  fairing  or  smoothing.  Also  shown  are  the  theoretical  boundaries 
as  calculated  by  the  method  of  characteristics.  Since  tolerances  in 
machining  the  test  nozzles  will  cause  errors  in  the  jet  initial  angle,  6  j, 
the  actual  initial  angle  is  shown  and  the  theoretical  boundaries  recalcu¬ 
lated  if  the  test  value  deviated  more  than  1  deg  from  the  desired  60  deg. 

The  experimental  data  generally  lie  slightly  outside  the  method-of- 
characteristics  solutions,  as  would  be  expected  for  a  viscid  fluid. 

However,  little  effect  of  gas  stagnation  temperature  or  pressure  is 
observed,  indicating  little  effect  of  nozzle  boundary  layer,  gas  con¬ 
densation,  or  changes  in  the  boundary  viscous  mixing.  The  exception 
is  the  data  for  N2  at  high  stagnation  pressures. 

The  boundary  increase  shown  for  the  high  stagnation  pressure  N2 
is  real,  as  previously  shown  by  the  pressure  probe  data  in  Fig.  6a,  and 
not  just  an  apparent  change  caused  by  the  flow  visualization  technique. 

Such  a  change  could  be  caused  by  gas  condensation  in  the  test  nozzles 
or  in  the  jet  plume.  The  very  large  rates  of  gas  expansion  in  the  noz¬ 
zles  and  in  the  plumes  would  give  a  substantial  amount  of  supersatura¬ 
tion  before  condensation  would  occur.  However,  as  shown  in  Fig.  12, 
even  equilibrium  condensation  in  the  test  nozzles  is  possible  only  for 
the  CO2  at  300°F  and  condensation  with  the  100°F,  100-psia  stagnation 
test  conditions  would  occur  far  before  condensation  at  the  1000°F, 

500-psia  conditions.  The  close  agreement  of  the  nitrogen  boundaries 
at  low  stagnation  pressure  test  conditions  and  for  all  conditions  with 
CO2  shows  that  gas  condensation  in  the  jet  plume  either  does  not  occur 
or  has  no  effect  on  the  measured  jet  boundaries. 

An  analysis  was  made  of  the  experimental  N2  data  for  possible  effects 
of  the  nozzle  boundary  layer  or  changes  in  the  jet  boundary  mixing  region. 
The  maximum  change  in  the  nozzle  exit  flow  properties  that  can  be  caused 
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by  the  displacement  thickness  of  a  turbulent  boundary  layer  will  cause 
only  a  0,  5-deg  increase  in  jet  boundary  initial  angle,  6j.  This  change 
is  not  sufficient  to  cause  the  boundary  location  increase  shown  for  the 
high  stagnation  pressures.  However,  an  analysis  of  the  deviation  of  the 
experimental  boundaries  from  the  method- of-characteristics  solutions 
shows  excellent  agreement  when  the  data  are  correlated  by  the  Reynolds 
number  at  the  nozzle  exit.  This  correlation  is  shown  in  Fig.  13,  which 
also  includes  the  data  with  He  and  CO2.  The  Reynolds  number  is  based 
on  the  length  along  the  nozzle  wall  from  the  throat  to  the  exit  and  the 
viscosity  data  of  Ref.  20  for  He  and  of  Ref.  21  (extrapolation  required 
for  lowest  temperatures)  for  N2  and  CO2.  The  critical  Reynolds  num¬ 
ber  indicated  in  Fig.  13  is  almost  identical  with  the  values  given  in 
Ref.  22  for  flat  plates  with  comparable  Mach  numbers.  A  somewhat 
higher  critical  Reynolds  number  would  be  expected  for  the  nozzle  flow 
because  of  the  favorable  pressure  gradient;  however,  the  choice  of  the 
proper  length  for  nozzle  flows  is  uncertain.  Also  shown  in  Fig.  13 
are  the  edges  of  the  turbulent  mixing  region  calculated  as  described  in 
Section  2.  3,  which  gives  good  agreement  with  the  data  at  the  higher 
Reynolds  numbers.  Although  it  would  be  attractive  to  attribute  the 
smaller  jet  boundaries  at  the  lower  Reynolds  numbers  to  a  laminar  jet 
boundary.  Ref.  23  would  indicate  laminar  boundary  thicknesses  com¬ 
parable  to  or  thicker  than  that  for  a  turbulent  boundary.  However  the 
velocity  profile  assumed  in  Ref.  23  may  not  be  valid  at  the  very  high 
Mach  numbers  of  the  experimental  jet  boundaries.  It  is  significant  that 
the  inviscid  method  of  characteristics  does  give  a  good  representation 
of  the  physical  jet  boundary,  although  insufficient  knowledge  of  the  vis¬ 
cous  mixing  process  and  the  factors  affecting  it  may  not  always  allow 
an  absolute  correction  for  the  mixing. 

In  general,  excellent  agreement' is  obtained  between  the  inviscid 
method-of-characteristics  solutions  and  the  experimental  data.  It  is 
therefore  believed  that  the  validity  of  this  program  with  the  inputs 
and  techniques  used  has  been  confirmed  and  that  a  comparison  of  simu¬ 
lation  parameters  can  now  be  made  using  the  method-of-characteristics 
calculations  for  the  jet  boundaries  rather  than  requiring  experimental 
measurements. 


4.3  ANALYSIS  OF  OTHER  SIMULATION  PARAMETERS 

Although  the  proposed  simulation  parameters  6j  and  yMi2/$i  do 
provide  good  simulation  for  most  cases  of  practical  interest,  they  do 
not^appear  to  be  the  completely  general  tool  that  one  desires.  There¬ 
fore  with  the  method-of-characteristics  program  confirmed  as  a 
method  of  calculating  the  physical  jet  boundary,  a  limited  analysis  of 
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some  of  the  other  proposed  simulation  parameters  was  made.  Among 
the  parameters  investigated  were: 

1.  yM j2/ /3j  and  6j 

2.  YMj2/j8j  and  Pj/Po,  (Ref.  3) 

3.  yM j 2 / |S j „  pj/p^andfij  (Ref.  1) 

Results  of  this  limited  analysis  are  shown  in  Fig.  14.  Of  the  sets 
of  parameters  studied  none  except  the  third  set  gave  results  as  good 
as  the  yM and  suggested  in  Ref.  4.  However  it  should  be  noted  . 
that  the  third  set  requires  the  matching  of  three  quantities  rather  than 
two,  and  is  therefore  a  more  restricted  simulation. 

The  use  of  the  third  set  of  simulation  parameters  shown  above 
implicitly  specifies  a  value  of  0]\j.  The  required  value  of  0jsf  is  there¬ 
fore  shown  in  Fig.  14c  which  presents  a  set  of  boundaries  calculated 
for  the  given  set  of  simulation  parameters. 


4.4  SIMULATION  PARAMETERS  S,  AND  Mi/y 

A  rapid  analysis  of  any  other  parameter  in  combination  with  6  j 
can  be  made  from  the  existing  method-of-characteristics  solutions  for 
6j  =  60  deg.  This  is  done  by  observing  the  variation  of  the  jet  bound¬ 
ary  radial  locations  at  a  fixed  axial  station  as  compared  to  changes  of 
the  parameter  under  study.  Ideally  then,  boundary  locations  for  all 
calculations  should  give  a  smooth  progression  for  the  changes  of  the 
selected  parameter.  A  summary  of  the  calculated  boundaries  from 
which  the  analysis  can  be  made  is  shown  in  Fig.  15. 

The  analysis  described  above  was  made  on  several  parameters 
involving  combinations  of  Mi  and  Y.  This  analysis  shows  excellent 
correlation  of  the  boundary  location  when  6j  and  Mi/y  are  used  as 
the  simulation  parameters.  A  slight  improvement  in  correlation  can 
be  gained  by  also  duplicating  0jq-  as  a  third  parameter.  Sample  results 
of  this  analysis  are  shown  in  Fig.  16.  Excellent  correlation  of  the 
boundary  locations  is  shown. 

Based  on  the  preceding  analysis,  additional  boundary  calculations 
were  made  by  the  method  of  characteristics  to  determine  the  range  of 
applicability  of  the  parameters  6j  and  Mi/y.  Calculations  were  made 
for  the  following  conditions: 
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1.  Jet  boundary  initial  angles,  6j,  of  45,  60,  and  75  deg. 

2.  Simulation  parameter  Mi/ y  values  of  5,  7,  10,  and 
13. 

3.  Specific  heat  ratio  of  1.  28  and  1.  667. 

4.  Nozzle  half-angle,  0]^,  of  10  and  20  deg. 

The  very  wide  range  of  conditions  covered  by  these  calculations  is 
shown  in  Table  II  and  in  Fig.  17. 

Results  of  the  calculations  are  shown  in  Figs.  18,  19,  and  20. 
Excellent  agreement  between  the  boundaries  is  obtained  over  the  very 
wide  range  of  conditions  investigated.  Some  improvement  is  obtained 
by  including  the  nozzle  half -angle,  @n,  as  a  third  simulation  parameter. 
Also  note  that  the  boundaries  are  presented  for  two  extreme  cases  of 
the  specific  heat  ratio,  1.  28  and  1.  667  only;  for  most  practical  applica¬ 
tions  using  diatomic  gases,  the  boundaries  would  show  even  better 
agreement. 


4.5  A  RAPID  METHOD  FOR  ESTIMATING  THE  INVISCID  BOUNDARY  LOCATION 

The  close  correlation  of  jet  boundaries  with  the  simulation  param¬ 
eters  6j,  M1/7,  and  0n  gives  rise  to  a  very  simple  method  of  rapidly 
estimating  the  inviscid  boundary  location.  Boundary  locations  at  fixed 
axial  stations  are  shown  in  Figs.  21  and  22  as  a  function  of  the  jet 
initial  angle,  6j,  and  the  parameter  Mj/t  for  nozzle  half-angles  of  10 
and  20  deg.  The  estimating  procedure  is  as  follows: 

1.  From  Eq.  (1)  the  initial  angle  of  the  jet,  6j,  is  calcu¬ 
lated.  This  establishes  the  slope  of  the  boundary  at 
the  nozzle  exit. 

2.  Calculate  M^/y  from  the  known  stagnation  pressure, 
ambient  pressure,  and  specific  heat  ratio. 

3.  The  boundary  location  r/re  can  then  be  located  at 
x/re  =  5,  10,  18,  and  25  from  Fig.  21  for  0n  =  10  deg 
or  Fig.  22  for  0jg  =  20  deg.  This  provides  up  to  five 
known  points  plus  the  initial  angle,  6j,  and  is  normally 
enough  for  an  accurate  boundary.  If  desired,  additional 
figures  similar  to  Figs.  21  and  22  can  be  constructed 
at  any  x/re  from  Figs.  18,  19,  and  20. 

4.  The  small  correction  for  other  values  of  0jj  can  be  made 
by  a  comparison  of  the  boundaries  at  10  and  20  deg. 
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SECTION  V 

CONCLUDING  REMARKS 


The  theoretical  and  experimental  investigation  of  jet  boundary 
simulation  parameters  for  underexpanded  jets  in  a  quiescent  atmos¬ 
phere  indicate  the  following  conclusions: 

1.  Experimental  jet  boundaries  are  in  excellent  agreement 
with  the  inviscid  jet  boundaries  as  determined  by  the 
method  of  characteristics. 

2.  Pindzola's  parameters  (6j  and  yM^/^)  give  good  boundary 
simulation  for  10  to  15  nozzle  exit  radii  for  most  applica¬ 
tions  using  air  (y  =  1.4)  to  simulate  a  hot  rocket  (y  =  1.  28). 

3.  However,  the  use  of  6j  and  Mi/y  provides  excellent  bound¬ 
ary  duplication  over  a  wider  range  of  conditions  and  to  much 
larger  axial  distances.  In  addition,  a  slight  improvement 

in  the  correlation  can  be  gained  by  also  duplicating  0n  as 
a  third  parameter. 

4.  A  procedure  developed  from  the  method- of- characteristics 
solutions  gives  a  rapid,  simple  method  for  estimating  the 
jet  boundary  shape  for  jets  with  boundary  initial  angles 
from  45  to  75  deg. 
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Fig.  1  18-lnch  Test  Cell  Installation 
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Fig.  4  Typical  Photographs  of  Experimental  Jet  Plumes  foryM1(//8i  **  14 
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Fig,  5  Typical  Jet  Density  Distribution 
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Fig.  6  Comparison  of  Boundaries  Determined  by  the  "Glow  Discharge"  Technique 
with  the  Pressure  Measured  by  a  Traversing  Probe 
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NOZZLE  SPHERICAL  AREA  RATIO,  A,/A 
a.  =  10 

Fig.  7  Range  of  Nozzle  Area  Ratio  and  Exit  Pressure  Ratio  for  Jets  with  Constant 
Values  of  yM*  /£,  and  S j  =  60  deg 
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Fig.  7  Continued 
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Fig.  8  Comparison  of  Jet  Boundaries  as  Determined  by  the  Latvala  Approximation 
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c.  =  18,  5j  -  60  deg 

Fig.  8  Continued 
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Fig.  9  Comparison  of  Jet  Boundaries  as  Determined  by  a 
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Fig.  10  Continued 
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Fig.  10  Continued 
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Fig.  10  Continued 
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Fig.  10  Continued 
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Fig.  11  A  Comparison  of  Experimentally  Determined  Jet  Boundaries  with  a  Method-of-Characteri sties 
Solution  for  --  18 
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Fig.  11  Continued 
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AXIAL  DISTANCE  TO  EXIT  RADIUS  RATIO,  x/r9 
d.  He,  0N  =  20  deg,  5j  -  60  deg 
Fig.  11  Concluded 
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Fig.  13  Effect  of  the  Nozzle  Exit  Reynolds  Number  on  the  Location  of  the  Experimental  Jet  Boundary  Relative  to 
the  Calculated  Inviscid  Boundary  at  x/re  =  20 
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Fig.  14  Comparison  of  Jet  Boundaries  for  Other  Proposed  Simulation  Parameters 
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Fig.  15  Concluded 
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BOUNDARY  MACH  NUMBER  TO  SPECIFIC  HEAT  RATIO,  M,/y 

a.  x/re  =  10 

Fig.  16  Correlation  of  Jet  Radii  from  Figure  15  at  a  Fixed  Axial  Distance  to  Exit  Radius  Ratio 
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b.  x/re  =  18 
Fig.  16  Concluded 
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a.  <5j  ==  75  deg 

Fig.  17  Range  of  Nozzle  Area  Ratio  and  Exit  Pressure  Ratio  for  Jets  with 
Constant  Values  af  M,/y 
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Fig.  18  Continued 
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Fig.  18  Continued 
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Fig.  18  Concluded 
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G.  M ,/y  =  13 

Fig,  19  Comparison  of  Jet  Boundaries  as  Determined  by  a  Method-of-Characteristics 
Solution  for  5j  =  60  deg 
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Fig.  19  Continued 
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Fig.  19  Continued 
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Fig.  19  Concluded 


68 


RADIAL  DISTANCE  TO  EXIT  RADIUS  RATIO,  r/r 


20 


a.  M \/y  =  10 

Fig.  20  Comparison  of  Jet  Boundaries  as  Determined  by  a  Method-of-Charocteri sties 
Solution  for  §j  -  45  deg 
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Fig.  20  Continued 
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Fig.  20  Concluded 
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Fig.  21  Effect  of  the  Boundary  Initial  Angle  on  the  Jet  Radius  at  a  Given  Axial  Distance 
for  Nozzles  with  =  10  deg 
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Fig.  21  Continued 
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Fig.  21  Continued 
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Fig.  21  Concluded 
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Fig.  22 


Effect  of  the  Boundary  Initial  Angle  on  the  Jet  Radius  at  a  Given  Axial  Distance 
for  Nozzles  with  0^  =  20  deg 
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Fig.  22  Continued 
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Fig.  22  Continued 


78 


RADIAL  DISTANCE  TO  EXIT  RADIUS  RATIO,  r/r, 


AEDC-TR-68-108 


40  50  60  70  80 

JET  BOUNDARY  INITIAL  ANGLE,  Sj,  deg 
d.  x/re  =  25 
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TABLE  1 

FLOW  PROPERTIES  OF  COMPUTED  JETS  MATCHING  THE 
SIMULATION  PARAMETERS  rtf /fa  AND  3j 


Case 

6j' 

deg 

/gl 

deg 

r 

M1 

Ag/A* 

“j 

Vp- 

1 

6 

0 

l 

0 

10 

1.28 

7.747 

8.677 

3.399 

345.5 

2 

20 

1.28 

7.747 

16.016 

3.920 

147.8 

3 

10 

1.40 

7.071 

2.887 

2.597 

222.4 

4 

| 

20 

1.40 

7.071 

4.585 

3.083 

106.1 

S 

10 

1.667 

5.912 

1.187 

1.566 

127.7 

6 

20 

1.667 

5.912 

1.533 

2.002 

68.2 

7 

l 

4 

10 

1.28 

10.891 

20.946 

4.152 

1806 

8 

20 

1.28 

10.891 

44.896 

4.836 

646.7 

g 

10 

1.40 

9.949 

4.786 

3.129 

921.9 

10 

20 

1.40 

9.949 

8.472 

3.740 

384.5 

n 

10 

1.667 

8.338 

1.441 

1.904 

396.0 

12 

20 

1.667 

8.338 

2.006 

2.406 

195.4 

13 

l 

8 

10 

1.28 

14.027 

38.124 

4.686 

7331 

14 

20 

1.28 

14.027 

90.935 

5.  509 

2304 

15 

10 

1.40 

12.818 

6.719 

3.489 

3010 

16 

20 

1.40 

12.818 

12 . 807 

4.201 

1143 

17 

10 

1.667 

10.751 

1.649 

2.114 

1003 

18 

20 

1.667 

10.751 

2.404 

2.672 

467.3 

19 

2 

2 

10 

1.28 

17.158 

58 . 573 

5.084 

24714 

20 

20 

1.28 

17.158 

151.41 

6.024 

7033 

21 

10 

1.40 

15.682 

8.549 

3.750 

8279 

22 

20 

1.40 

15.682 

17.197 

4.545 

2927 

23 

10 

1.667 

13.159 

1.816 

2.259 

2200 

24 

f 

20 

1.667 

13.159 

2.733 

2.861 

982.5 

80 


TABLE  II 

- 

FLOW  PROPERTIES  OF  COMPUTED  JETS  MATCHING  THE 

SIMULATION  PARAMETERS  M,/y  AND  S| 

Case 

6j* 

M  X/Y 

eN> 

Y 

M1 

VA* 

M. 

deg 

deg 

j 

25 

1 

r5 

5 

10 

1.28 

6.400 

2.728 

2.405 

26 

20 

1.28 

6.400 

4.135 

2.771 

27 

10 

1.667 

8.335 

1.047 

1.270 

28 

20 

1.667 

8.335 

1.265 

1.685 

29 

7 

10 

1.28 

8.960 

5.504 

3.015 

30 

20 

1.28 

8. 960 

9.418 

3.468 

31 

10 

1.667 

11.669 

1.152 

1.506 

32 

1 

20 

1.667 

11.669 

1.469 

1.934 

33 

1( 

3 

10 

1.28 

12 . 80Q 

11.123 

3.609 

34 

• 

20 

1.28 

12.800 

21.404 

4.171 

35 

10 

1.667 

16.670 

1.265 

1.686 

36 

20 

1.667 

16.670 

1.679 

2.141 

37 

IS 

I 

10 

1.28 

16.640 

17.541 

3.998 

38 

20 

1.28 

16.640 

36.465 

4.645 

39 

10 

1.667 

21  671 

1.341 

1.787 

40 

20 

1.667 

21.671 

1.820 

2.262 

41 

61 

0 

c 

K. 

> 

10 

1.28 

6.400 

5.230 

2.  972 

42 

20 

1.28 

6.400 

8.810 

3.412 

43 

10 

1.667 

8.335 

1.410 

1 . 940 

44 

20 

1.667 

8.335 

2.005 

2.406 

Vp- 


406.6 

217.3 

976.9 

541.5 

2201.0 

1029.0 

3580.0 

1933.0 

17527.0 

7144.0 

15806.0 

8227.0 

95410.0 

35450.0 

50404.0 

25621.0 

154.1 
73.2 

395.6 

195.1 
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oo 

to 


Case  M^/y 

deg 


45  60 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57  45 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 


TABLE  II  (Concluded) 


°NJ 

deg 

Y 

M1 

VA* 

"j 

Pi/Pa. 

J  “ 

10 

1.28 

8.960 

12.729 

3.723 

677.6 

20 

1.28 

8. 960 

25.054 

4.310 

268.8 

10 

1.667 

11.669 

1.717 

2.175 

1373.0 

20 

1.667 

11.669 

2.536 

2.751 

628.4 

10 

1.28 

12.800 

30.918 

4.497 

4345.0 

20 

1.28 

12.800 

70.991 

5.268 

1431.0 

10 

1.667 

16.670 

2.005 

2.406 

5694.0 

20 

1.667 

16.670 

3.116 

3.057 

2434.0 

10 

1.28 

16 . 640 

55.084 

5.026 

20419.0 

20 

1.28 

16.640 

140.680 

5.  948 

5898.0 

10 

1.667 

21.671 

2.202 

2.543 

17450.0 

20 

1.667 

21.671 

3.536 

3.250 

7101.0 

10 

1.28 

6.400 

11.930 

3.669 

47.7 

20 

1.28 

6.400 

23 . 200 

4.243 

19.4 

10 

1.667 

8.335 

2.462 

2.707 

13.0 

20 

1.667 

8.335 

4.085 

3.475 

50.6 

10 

1.28 

8. 960 

36.667 

4.650 

160.9 

20 

1.28 

8. 960 

86.839 

5.463 

51.0 

10 

1.667 

11.669 

3.224 

3.109 

399.5 

20 

1.667 

11.669 

5.844 

4.061 

136.0 

10 

1.28 

12.800 

113.81 

5.732 

765.  9 

20 

1.28 

12.800 

335.63 

6.887 

185.0 

10 

1.667 

16.670 

4.084 

3.474 

1479.0 

20 

1.667 

16.670 

8.035 

4.625 

443.7 
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